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Abstract 

Introduced  is  a  lattice-gas  with  long-range  2-body  interactions.  An 
effective  inter-particle  force  is  mediated  by  momentum  exchanges.  There 
exists  the  possibility  of  having  both  attractive  and  repulsive  interactions 
using  finite  impact  parameter  collisions.  There  also  exists  an  interesting 
possibility  of  coupling  these  long-range  interactions  to  a  heat  bath.  A 
fixed  temperature  heat  bath  induces  a  permanent  net  attractive  interpar¬ 
ticle  potential,  but  at  the  expense  of  reversibility.  Thus  the  long-range 
dynamics  is  a  kind  of  a  Monte  Carlo  Kawasaki  updating  scheme.  The 
model  has  a  PpT  equation  of  state.  Presented  are  analytical  and  numeri¬ 
cal  results  for  a  lattice-gas  fluid  governed  by  a  nonideal  equation  of  state. 
The  model’s  complexity  is  not  much  beyond  that  of  the  FHP  lattice-gas. 
It  is  suitable  for  massively  parallel  processing  and  may  be  used  to  study 
critical  phenomena  in  large  systems. 


1  Introduction 

Nonideal  fluids,  with  dynamics  governed  by  reversible  physical  laws,  undergo 
phase  transitions.  This  fact  about  fluids  indicates  the  possibility  that  lattice-gas 
fluid  models,  with  dynamics  governed  by  reversible  rules  [1],  may  also  undergo 
phase  transitions.  The  Ising  model  is  the  most  well  known  computational  model 
with  an  order-disorder  transition.  Reversible  Ising  models  using  energy  bankers, 
in  a  microcanonical  ensemble,  are  known  [2,  3],  Yet  it  is  an  open  question  as  to 
whether  or  not  there  exists  a  reversible  momentum-conserving  lattice-gas  model 
of  a  multiphase  fluid. 
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In  molecular  dynamics  one  simulates  a  many-body  system  of  particles  with 
continuous  interaction  potentials  where  the  particles  have  continuous  positions 
and  momenta.  In  lattice-gas  dynamics  the  particles’  positions  and  momenta  are 
discrete  and  motion  is  constrained  to  a  spacetime  lattice.  Interparticle  potentials 
can  be  modeled  by  including  long-range  interactions  in  the  lattice-gas  dynamics 
with  a  discrete  momentum  exchange  between  particles.  The  use  of  momentum 
exchange  was  introduced  by  Kadanoff  and  Swift  in  a  Master-equation  approach 
[4].  The  use  of  negative  momentum  exchanges  in  long-range  interactions  was 
first  done  in  a  lattice-gas  model  by  Appert  and  Zaleski  [5].  This  nonthermal 
model  has  a  liquid-gas  coexistence  phase;  there  is  a  Pp  equation  of  state.  A 
method  for  modeling  interparticle  potentials  using  only  local  interactions  was 
introduced  by  Chen  et  al.  [6].  There  is  an  Ising  interaction  between  the  rest 
particles  of  the  FHP  lattice-gas  model[7]  with  rest  particles  at  the  neighbor¬ 
ing  sites.  In  this  way  a  local  configurational  energy  is  associated  with  the  rest 
particles.  Speed  one  particles  can  transition  to  a  rest  state  with  a  certain  Boltz¬ 
mann  probability,  e_/3AE.  The  inverse  transition  is  also  possible  and  the  model 
obeys  detailed-balance.  Chen  et  al.  observed  an  order-disorder  transition  as 
the  system  has  a  nonideal  equation  of  state.  It  is  a  purely  local  model  being  a 
combination  of  a  lattice-gas  automaton  and  a  Monte  Carlo  Ising  lattice  gas.  In 
this  way  momentum  conservation  is  added  to  Ising  dynamics  so  the  model  can 
therefore  be  used  to  view  the  kinetics,  even  near  the  critical  point. 

The  lattice-gas  model  with  long-range  interactions  presented  here  is  a  finite 
temperature  extension  of  Appert  and  Zaleski ’s  zero  temperature  model.  The 
first  ingredient  added  is  repulsive  long-range  collisions.  Both  discrete  negative 
and  positive  momentum  exchanges  occur  between  particles.  The  second  ingre¬ 
dient  added  is  a  finite  temperature  heat-bath,  that  is  a  heat-bath  with  a  certain 
non-zero  fractional  occupation.  It  is  possible  to  bias  the  finite  impact  parame¬ 
ter  collisions  so  there  is  a  net  attractive  interparticle  potential.  This  is  done  by 
coupling  the  long-range  collisions  to  a  heat  bath  —  attractive  collisions  cause 
a  transition  from  a  high  potential  energy  state  to  a  low  one  and  emit  units  of 
heat  whereas  repulsive  collisions  cause  the  opposite  transition  and  absorb  heat. 
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When  a  disordered  lattice-gas  state  is  in  contact  with  a  low  temperature  bath, 
spinoidal  decomposition  occurs.  The  fractional  occupation  of  the  heat-bath  de¬ 
termines  the  likelihood  of  long-range  interactions  in  a  simple  way.  The  phase 
separation  occurs  if  the  heat  bath  filling  fraction  is  held  fixed1.  Its  nonideal 
PpT  equation  of  state  is  derived  and  compared  to  numerical  simulation. 

When  the  lattice-gas  dynamics  is  strictly  reversible,  there  exists  an  inherent 
limitation  that  the  phase  separation  process  can  occur  only  for  a  short  period 
of  time.  The  lattice-gas  fluid  quickly  becomes  a  neutral  fluid,  with  finite  impact 
parameter  collisions.  So  added  to  the  usual  FHP  type  on-site  collisions  are  an 
equivalent  set  of  finite  impact  parameter  collisions.  Balancing  the  interactions 
ensures  detailed  balance,  and  in  the  context  of  the  multiphase  model  presented 
below,  this  is  like  the  infinite  temperature  limit.  The  appendix  contains  a 
description  of  the  reversible  lattice-gas  with  balanced  attractive  and  repulsive 
collisions  and  a  numerical  result  illustrating  the  characteristic  transient  time. 

This  paper  is  organized  into  three  main  sections.  §2  very  briefly  describes 
the  local  particle  dynamics  of  the  lattice-gas  method.  §3  describes  the  long- 
range  lattice-gas  and  offers  a  simple  theoretical  result  in  the  Boltzmann  limit. 
Finally,  §4  presents  some  numerical  results  obtained  with  the  model.  A  closing 
discussion  of  the  main  points  of  this  paper  is  given  in  §5.  The  appendix  contains 
a  formal  construction  of  a  long-range  lattice-gas  with  a  single  species  of  particles 
obeying  detailed  balance. 

2  Lattice- Gas  Automata 

An  extremely  abridged  description  of  local  lattice-gas  dynamics  is  given  here 
since  descriptions  can  be  found  elsewhere  [8,  9],  Particles,  with  mass  m,  prop¬ 
agate  on  a  spacetime  lattice  with  N  spatial  sites,  unit  cell  size  l,  time  unit  r, 
with  speed  c  =  ■£.  A  particle’s  state  is  completely  specified  at  some  time,  t, 
by  specifying  its  position  on  the  lattice,  x,  and  its  momentum,  p  =  meea  with 
lattice  vectors  ea  for  a  =  1,2, ... ,2?.  The  particles  obey  Pauli  exclusion  since 

'Note  that  holding  the  heat  bath  at  fixed  fractional  occupation  below  0.5  breaks  detailed 
balance. 


3 


only  one  particle  can  occupy  a  single  state  at  a  time.  The  total  number  of  con¬ 
figurations  per  site  is  2s.  The  total  number  of  states  available  in  the  system  is 
2BN .  The  lattice-gas  cellular  automaton  equation  of  motion  is 

na(x  +  /ea,  t  +  r)  =  na(x,t)  +  Qa(n(x,t)),  (1) 

where  the  particle  occupation  variable  and  collision  operator  are  denoted  by  na 
and  ft a,  respectively. 

For  a  two-dimensional  hexagonal  lattice,  the  spatial  coordinates  of  the  lat¬ 
tice  sites  may  be  expressed  as  follows  Xij  =  {^3>  *  —  5  mod  2j\  where  i  and  j 
are  rectilinear  indices  that  specify  the  memory  array  locations  used  to  store  the 
lattice-gas  site  data.  Given  a  particle  at  site  (i,j),  it  may  be  shifted  along  vector 
r  =  reB  to  a  remote  site  ( i',j')a  by  the  following  mapping:  (i  +  —  mod2l?  mod  2r,j  T  r)x  , 

(i  -  5  -  mod2j  mod  2r,j  T  r)2  5,  (i  T  r,j)3  6.  These  streaming  relations  are 
equivalent  to  memory  address  offsets.  The  modulus  operator  is  base  2  because 
even  and  odd  rows  must  be  shifted  as  a  hexagonal  lattice  is  embedded  into  a 
square  lattice. 

3  Long-Range  2-Body  Interactions 

An  interparticle  potential,  V(x  —  x'),  acts  on  particles  spatially  separated  by 
a  fixed  distance,  x  —  x'  —  2r.  An  effective  interparticle  force  is  caused  by  a 
non-local  exchange  of  momentum.  Momentum  conservation  is  violated  locally, 
yet  it  is  exactly  conserved  in  the  global  dynamics. 

For  the  case  of  an  attractive  interaction,  there  exists  a  bound  states  in  which 
two  particles  orbit  one  another.  Since  the  particle  dynamics  are  constrained 
by  a  crystallographic  lattice  we  expect  polygonal  orbits.  In  figure  la  we  have 
depicted  two  such  orbits  for  a  hexagonal  lattice-gas.  The  radius  of  the  orbit  is  r. 

Two-body  finite  impact  parameter  collisions  are  depicted  in  figures  lb  and  lc. 

Momentum  exchanges  occur  along  the  principle  directions.  The  time-reversed 
partners  of  the  collisions  in  figures  lb  are  included  in  the  model.  The  interaction 
potential  is  not  spherically  symmetric,  but  has  an  angular  anisotropy.  In  general, 
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Figure  1:  (a)  Bound  state  orbits  where  the  dotted  path  indicates  the  particle’s  closed 

trajectory;  (b)  |Ap|  =  1  with  one  unit  of  angular  momentum,  a  counter-clockwise  attractive 
collision  and  its  repulsive  conjugate;  and  (c)  |Ap|  =  2  with  zero  angular  momentum  collision 
conjugates.  Not  included  in  the  figure  are  the  time-reversed  partners  of  (b). 
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it  acts  only  on  a  finite  number  of  points  on  a  shell  of  radius  r.  The  number  of 
lattice  partitions  necessary  per  site  is  half  the  lattice  coordination  number,  since 
two  particles  lie  on  a  line.  Though  microscopically  the  potential  is  anisotropic, 
in  the  continuum  limit  numerical  simulation  done  by  Appert,  Rothman,  and 
Zaleski  indicates  isotropy  is  recovered  [10]. 

Constraint  equations2  for  momentum  conservation  and  parallel  and  perpen¬ 
dicular  momentum  exchange  are  respectively 

ea  -  ep  +  ep  -  e„  =  0  (2) 

(eQ  -  e/3  -  +  e„)  •  r  =  2A p  (3) 

(eQ  -  ep  -eM  +  e„)  x  r  =  0  (4) 

where  A p  is  the  momentum  change  per  site  due  to  long-range  collisions.  The 
sum  and  difference  of  (2)  and  (3)  reduce  to 

(®a)j  —  (®/3)y  =  (®p)y  —  (ev)y  =  —  Ap.  (5) 

The  possible  non-zero  values  of  a  site’s  momentum  change  may  be  A p  —  ±1 
and  ±2.  As  mentioned,  the  cases  for  A p  <  0  led  to  bound  states  with  angular 
momentum  0  and  1.  To  satisfy  (5),  consider  the  case  where  (e„)„  =  —  (eM)j, 
and  (e(j)y  —  —  (e„)y.  3  The  possible  collisions  where  f  —  e 3  are  depicted  in 
figure  1.  The  reversible  interactions  are  2-body  collisions  with  a  finite  impact 
parameter  of  2r.  For  r  =  0,  they  reduce  to  the  2-body  collisions  in  the  FHP 
lattice-gas:  the  |Ap|  =  1  collisions  reduce  to  rotations  of  momenta  states, 
and  the  |Ap|  =  2  collisions  reduce  to  the  identity  operation. 

Let  V[d(a:)]  represent  the  potential  energy  due  to  long-range  interactions, 
where  d(x)  is  the  probability  of  finding  a  particle  at  position  x.  For  a  2-body 
interaction  to  occur  at  x  and  x' ,  one  must  count  the  chance  of  having  two 
particles  and  two  holes  at  the  right  locations,  so  in  the  Boltzmann  limit  one 

2We  are  simplifying  this  development  by  assuming  a  single  speed  lattice-gas.  Consequently 
we  do  not  have  to  explicitly  write  a  term  to  conserve  energy  since  here  energy  conservation 
follows  by  default. 

Alternatively  one  could  have  chosen  (eQ)y  =  (e„)„  and  (ep)y  =  (e^)j,. 
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Figure  2:  Potential  energy  versus  particle  density  in  the  zero  temperature  limit.  Letting 
d  — ►  V  —  f  i  then  V(ip)  =  —  ^ip2  +  ifi4  which  has  a  Landau-Ginsburg  form. 

may  write  a  probability  of  collision,  P(x,x';v,v'),  as 

P(x,x’;v,v')  =  d(x)(  1  —  d(a;))d(x')(l  —  d(x'))J(|w'  —  v|)  (6) 

and  if  the  system  is  uniformly  filled,  this  simplifies  to 

P(x,  x'\  v,  v')  =  d2(  1  -  d)2SV'V.  (7) 

Letting  m,  c,  r  =  |x  —  x'|,  and  l  denote  the  particle  mass,  particle  velocity,  2- 
body  interaction  range,  and  lattice  cell  size,  one  may  write  the  potential  energy 
as 

V(d)  =  ^m<?B  (y)  P{x,x'\v,v')  (8) 

=  ^ mc2B  (y)  d2(  1  -  d)2SV’V,  (9) 

where  the  value  of  the  coefficient  a  depends  on  the  magnitudes  of  momenta 
exchanged.  Here  d  ranges  from  0  to  1,  and  is  just  the  particle  filling  fraction. 
V(d)  has  two  minima,  at  d  =  0  and  d  =  1,  see  figure  2. 
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Figure  3:  Examples  of  long-range  collisions  that  locally  conserve  mass,  momentum,  and 
energy.  |Ap|  =  ±1  interactions  along  the  ro-direction  coupled  to  a  heat  bath:  (a)  attractive 
case;  and  (b)  its  adjoint,  repulsive  case.  Transitions  probabilities  for  these  collisions  have  a 
Monte  Carlo  form  as  they  are  biased  by  the  density  of  heat-bath  particles,  rendered  here  with 
wavy  lines.  Head  of  the  gray  arrows  indicates  particles  entering  the  sites  at  1*0  and  —  ro  at 
time  t.  Tail  of  the  black  arrows  indicates  particles  leaving  those  sites  at  time  t  +  T. 


One  may  consider  a  slightly  more  complicated  interaction,  where  the  2-body 
collisions  are  coupled  to  a  second  kind  of  particle  whose  filling  fraction  is  denoted 
by  h.  In  the  slightly  more  complicated  interaction,  the  form  of  V  (d)  given  still 
holds,  but  only  for  h  =  0.  Here  is  the  slightly  more  complicated  version  of 
things.  The  complete  form  of  the  interaction  energy  is 

V(d,  h)  =  |mc2S  [d2(l  -  d)2(l  -  h)2  -  d2(l  -  dfh2)  .  (10) 

The  first  term  is  two  d’s  transitioning  to  a  lower  configurational  energy  state 
and  thus  emitting  two  h’s  to  conserve  energy.  The  second  term  is  two  d’s 
transitioning  to  a  higher  configurational  energy  state  by  absorbing  two  h’s. 
Local  conservation  of  momentum  and  energy  is  recovered.  For  convenience  we 
write 

V(d,  h )  =  ^mc2B  [d2(l  -  d)2(l  -  2h)}  .  (11) 

From  (11),  V(d,h)  =  0  for  h  =  A.  Since  the  hath-bath  particles  are  fermi-dirac 
distributed,  we  define  the  effective  temperature  as  kgT  =  e0(log  ^jj®)-1,  and 
h  =  |  corresponds  to  T  =  oo  and  h  =  0  corresponds  to  T  =  0.  Numerical 
simulation  corroborates  this.  The  pressure,  p,  in  the  gas  is  written 

p{d,h)=mc2Bd  +  V{d,h)  (12) 
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where  cs  is  the  sound  speed.  This  is  the  non-ideal  equation  of  state  that  is 
responsible  for  the  liquid-gas  phases  observed  in  numerical  simulations  of  this 
system. 

4  Simulation  Results 

The  dynamics  of  the  model  in  contact  with  a  fixed  temperature  heat  bath  is 
tested  by  numerical  simulation.  A  coarse-grained  mass  frequency  distribution  is 
measured  after  the  system  has  evolved  for  a  fixed  amount  of  time.  The  lattice- 
gas  is  initialized  with  a  random  configuration  and  allowed  to  evolve  for  500  time 
steps  for  several  bath  filling  fractions:  30%  to  20%,  19%,  18%,  5%,  and  0%. 
Resulting  system  snapshots  are  illustrated  in  figure  4. 

If  the  lattice-gas  is  above  the  transition  temperature,  the  particles  are  uni¬ 
formly  spread  over  the  lattice.  As  the  system  evolves  while  in  contact  with  a 
finite  temperature  heat  bath  coarse-grained  4x4  block  averages  over  a  256  x  256 
lattice  are  taken  over  the  lattice-gas  number  variables  to  produce  a  mass  fre¬ 
quency  distribution  for  a  large  number  of  temperatures  and  the  liquid  and  gas 
densities  are  found.  A  mass  frequency  distribution  obtained  by  this  coarse¬ 
grained  block  averaging  procedure  is  a  Gaussian  with  its  mean  located  exactly 
at  the  initial  particle  density.  A  normalized  Gaussian  fit,  centered  at  particle 
density  0.3,  is  shown  in  figure  5  as  is  cross-section  plots  of  the  distribution  at 
different  temperatures.  As  the  temperature  decreases,  the  distribution  widens 
and  becomes  bimodal.  The  mean  of  the  low  density  peak  gives  the  gas  phase 
density  and  the  mean  of  the  high  density  peak  gives  the  liquid  phase  density. 
The  order  parameter  for  the  liquid-gas  transition  is  the  difference  of  the  liquid 
and  gas  densities,  ip  =  pi  —  pG  [11]- 

If  the  heat-bath  temperature  is  held  constant,  the  dynamics  is  no  longer 
reversible.  The  ordered  phase  persists  and  the  simulation  method  becomes 
like  a  Monte  Carlo  Kawasaki  updating  scheme  (i.e.  the  exchange  of  randomly 
chosen  spins).  However,  using  this  long-range  interaction  method,  momentum 
is  exactly  conserved  and  kinetic  information  retained.  Therefore,  the  dynamical 
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Figure  4:  Several  lattice-gas  configurations  obtained  after  500  iterations  starting  from  ran¬ 
dom  initial  256  x  256  configurations  with  30%  particle  filling.  The  six  configuration  are  coupled 
to  a  heat-bath  with  0,  5,  18,  19,  20,  and  30%  heat  bath  filling. 


Density 


Figure  6:  Liquid-gas  coexistence  curve  determined  from  the  particle  mass  frequency  distri¬ 
bution  for  different  heat  bath  temperatures. 
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evolution  of  the  finite  temperature  multiphase  system  is  accessible,  even  near  the 
critical  temperature.  Figure  6  shows  a  comparison  of  numerical  simulation  data 
obtained  by  this  procedure  to  an  analytical  calculation  done  in  the  Boltzmann 
limit  by  analytically  carrying  out  a  Maxwell  construction.  The  Gibbs  free  energy 
of  the  lattice-gas  can  be  written  analytically  since  the  pressure’s  dependence  on 
density  and  temperature  is  known.  With  (12)  and  defining  the  free  energy  as 

(13> 

a  Maxwell  construction  then  correctly  predicts  the  liquid  and  gas  phase  densi¬ 
ties  at  any  given  heat-bath  temperature,  see  figure  6.  This  mean  field  type  of 
calculation  itself  is  very  interesting  and  is  a  good  example  of  how  analytical  cal¬ 
culations  are  possible  in  simple  discrete  physical  models  like  lattice-gases.  The 
result  shown  in  figure  6  is  similar  to  the  order  parameter  curve,  magnetization 
versus  temperature,  of  an  Ising  model.  Spin  up,  (M)  =  +1,  and  spin  down, 
(M)  =  —  1,  domains  are  analogous  to  the  liquid  and  gas  phases. 

5  Discussion 

The  model  is  a  simple  discretization  of  molecular  dynamics  with  interparticle 
potentials.  Because  of  the  model’s  small  local  memory  requirement,  the  dynam¬ 
ics  of  large  systems  can  be  implemented  on  a  parallel  architecture,  as  has  been 
done  on  the  cellular  automata  machine  CAM-8  [12]. 

The  main  points  of  this  paper  are: 

1.  Coupling  the  particle  dynamics  to  a  fixed  temperature  heat-bath  sets  the 
transition  probabilities  and  causes  a  net  attractive  interparticle  potential  in  the 
macroscopic  limit.  The  heat-bath  is  comprised  of  a  set  of  lattice-gas  particles 
encoding  a  unit  of  heat.  The  heat  bath  density,  h,  controls  the  heat-bath’s 
temperature  by  the  fermi-dirac  distribution,  ksT  =  e0(log  ijp)-1.  The  system 
possesses  a  nonideal  PpT  equation  of  state.  With  the  model  in  contact  with  a 
fixed  heat-bath,  one  should  classify  it  as  finite  temperature  model  with  heat-bath 
dynamics  similar  to  a  Monte  Carlo  Kawasaki  updating,  yet  retaining  essential 
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kinetic  features.  It  is  a  step  toward  a  more  complete  long-range  lattice-gas  that 
preserves  an  interaction  energy. 

2.  The  equation  of  state  is  known  in  the  Boltzmann  limit.  A  Maxwell 
construction  predicts  the  liquid-gas  coexistence  curve.  Van  der  Waal  coefficients 
could  be  determined  to  map  the  simulation  on  to  a  particular  physical  liquid-gas 
system. 

3.  With  the  single  speed  particles  species  coupled  to  a  heat-bath  in  the  fash¬ 
ion  described  above,  there  is  an  imperfect  tracking  of  interaction  energy.  This  is 
a  limitation  of  this  reversible  lattice-gas.  An  improved  version  of  this  reversible 
lattice-gas  model  with  long-range  interactions  perhaps  could  be  implemented 
(likely  with  a  species  of  “bound  state”  particles  included)  so  that  simulations 
are  carried  out  in  a  microcanonical  ensemble  analogous  to  Ising  models  with 
auxiliary  demons  introduced  by  Creutz  [3],  and  Toffoli  and  Margolus  [2]. 

It  is  hoped  that  a  lattice-gas  model,  of  the  kind  presented  here,  will  become  a 
valuable  new  tool  for  analytically  and  numerically  studying  the  dynamic  critical 
behavior  of  multiphase  systems. 
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A  Reversibility-Neutrality  Statement 

The  main  thrust  of  this  paper  has  been  to  describe  a  momentum  conserving 
multiphase  lattice-gas.  In  this  appendix  a  reversible,  momentum-conserving 
lattice-gas  with  long-range  interactions  is  analysized.  A  formalism  is  presented 
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for  describing  the  unitary  evolution  of  a  lattice-gas  with  long-range  interactions. 

For  any  reversible  computational  model  a  unitary  operator  maps  the  compu¬ 
tational  state  at  some  time  to  the  state  at  the  next  time  iteration.  This  unitary 
matrix  can  be  expressed  as  the  exponential  of  a  hermitian  operator,  a  kind  of 
computational  Hamiltonian4. 

Using  the  notation  of  multiparticle  quantum  mechanical  systems  in  the  sec¬ 
ond  quantized  number  representation  [14],  all  states  of  the  system  are  enumer¬ 
ated  sequentially  by  ax  =  1 . . .  ArtotM  where  Ntol^  —  BN.  For  each  a  there  is 
an  associated  site  x  that  is  indicated  by  a  subscript.  Denote  the  vacuum  state 
of  the  system  by  |0)  where  all  n0l  =  0  for  all  states  a.  Using  creation  and  an¬ 
nihilation  operators  a^ai  and  aax  to  respectively  create  and  destroy  a  particle 
in  state  a  at  lattice  position  x,  any  arbitrary  system  configuration  tp  with  P 
particles  can  be  formed  by  their  successive  application  on  the  vacuum 

w=nat“,i°>’  (i4> 

?= i 

where  particle  one  is  in  state  «i  at  lattice  node  Xi,  particle  two  is  in  state 
at  lattice  node  X2,  etc.  The  number  operator  is  nQi  =  a* aiaax.  To  completely 
specify  the  dynamics  the  local  anticommutation  relations  are  required.  Since 
there  is  exclusion  of  boolean  particles  at  a  single  momentum  state,  we  have 


{aa,a^a}  =  1  (15) 

{«q,S«}  =  0  (16) 

{attt,ata}  =  0.  (17) 

However,  the  boolean  particles  are  completely  independent  at  different  momen¬ 
tum  states,  and  so  the  nonlocal  operators  commute 

=  o  (is) 

[aa,ap]  =  0  (19) 

[a^aV]  =  o  (20) 


4 This  mathematical  construction  is  similar  to  a  quantum  mechanical  description  [13]. 
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for  a  t^/3. 

A  unitary  evolution  operator  that  describes  the  complete  evolution  of  the 
lattice- gas  may  be  partitioned  into  a  streaming  and  collisional  part,  U°  and 
tjlnt  respectively.  The  full  system  evolution  operator  is  a  product  of  these  two 
operators  U  =  U°U  nt.  The  operator  U°  is  constructed  using  a  unitary  single 
exchange  operator  denoted  by  X'aP'  ■  All  permutations  of  single  boolean  par¬ 
ticle  states  may  be  implemented  by  successive  application  of  this  momentum- 
exchanger.  We  will  use  the  same  symbol,  x^/sm  to  denote  the  permutations 
between  state  a  at  site  x  and  states  (3 '  at  site  x'.  We  wish  to  construct  x^ 
from  the  boolean  lattice-gas  creation  and  annihilation  operators. 

We  require  that  x^  is  unitary,  (x^)2  =  1,  that  it  conserve  the  number  of 
particles,  [x(1*,  N]  =  0,  and  that  x^|  0)  =  |  0).  It  has  the  form 

Xa^  =  "t"  0^(3'  d”  1  ad  ad  0>  0,  o}  p,  dp' daQ^  (21) 


This  can  be  written  in  the  form 


(22) 


For  a  set  of  N  states,  {1,2,...,  TV  —  1,  TV},  two  Cn  rotation  operators  can  be 
implemented 

N-l 


<;(!)  *(« 
Xn-1,nXN-2,N-1 

..*(1)  -  TT  *(1) 

Xl2  —  11  XN~i,N-i+l 
i=l 

(23) 

X12  A23  Xn-i,n 

AT-1 

-  n  *is» 

*= 1 

(24) 

that  are  rotations  by  ±555°.  Suppose  we  pick  a  subspace  to  be  the  set  of  states, 
Pa,  with  momentum  mcea.  Then  following  our  construction,  we  have  found  a 
method  to  implement  a  unitary  streaming  operator,  Sa,  along  direction-a 


Sa 


IT  ■ 


The  free  part  of  the  evolution  is  then  simply 

B  B 


0  =  n  s. = n  n  «Sw 

a=l  o=l  {a,0}eVa 


(25) 


(26) 
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The  corresponding  kinetic  energy  part  of  the  Hamiltonian  to  leading  order  is 

H  =  'y  '  o}aiapx,  +  •  (27) 

(**'> 

where  the  sum  is  over  all  bonds  of  the  lattice  taken  over  partitions  along  principle 
lattice  directions  and  where  for  brevitity  the  following  short-hand  notation  is 

used:  £f=i  T,{a,0)eVa  -*■  £(**')  when  *  +  lea  -*■  *'• 

The  operator  U*  '  is  constructed  using  a  unitary  double  exchange  operator, 
denoted  by  X^n’u'  >  a  generalization  of  the  single  exchange  operator.  All  permu¬ 
tations  of  two  boolean  particles  may  be  implemented  by  successive  application 
of  this  momentum-exchanger,  where  the  permutations  for  particle  one  occurs 
between  states  a  and  0  at  site  x  and  for  particle  two  between  states  y!  and  i/ 
at  site  x'. 

We  require  that  is  unitary,  (y^)2  =  1,  that  it  conserve  the  number 
of  particles,  [x(2\fV]  =  0,  and  that  0)  =  |  0).  A  relation  identical  to  (22) 
exists  for  the  double  boolean  exchange  operator 


(2)  ,  o<T<2)  , 

\L'„'  =  1  -  2N.  =  e  **'=»»' . 


(28) 


Let  us  assume  we  have  a  two-particle  state  |'ycr/)  =  a)lxc0 Vx,\  0),  where 
7  =  a  or  0,  and  o'  =  y!  or  v' .  It  has  the  form 

X(x0nf  i/r  =  (0^CL0(0  00(^0^  1  (0 00,0(0 i/t  Oy*  0,^(0  ccO'^l'  (0  (0  ^0,^(0  ^>0,^0, 0(0  0Q,i/'Q^  , 

(29) 

Suppose  we  pick  a  subspace  to  be  the  set  of  states,  Va,  where  moment  exchanges 
±mceQ  can  occur  between  two  particle  pairs.  Then  a  unitary  collision  opera¬ 
tor,  CQ,  in  this  subspace  (with  momentum  exchanges  along  a  principle  lattice 
direction)  is 

{a,0  n,v}€Va 

The  interaction  part  of  the  evolution  is  then  simply 

B  B 

c'  -  n = n  n  <3» 

a=l  a=l  {a,(3  n,v}€Va 
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The  corresponding  potential  energy  part  of  the  Hamiltonian  to  leading  order  is 


H  —  ^  ^ axQ'/3xQ^ nxtQ'l/xi  ^ f3x^ax^ uxtafjixi  "T  *  *  *  j  (32) 

{xx,\vvt) 

B 

where  for  brevity  the  following  short-hand  notation  is  used:  ^a=i  fi  ^  „}<=y  — > 

Y2(xx‘-vv')  when  x  +  rea  — »•  x  and  x  —  rea  — >  x' .  (27)  and  (32),  imply  to  leading 
order,  the  full  lattice-gas  Hamiltonian  will  have  the  form 

H  =  H°+Hln,-I - -  ^2  a) axapxl  +  £  &Kx&px&^x,&ux,  +c.c.  +  -  •  •  (33) 

(xxf)  (ll'jw') 

The  complex  conjugate  terms  arise  in  (33)  because  of  reversibility  and  ensure 
the  hermiticity  of  the  Hamiltonian.  If  the  interaction  term  in  the  Hamiltonian 
covers  all  possible  attractive  interactions,  its  complex  conjugate  then  covers  all 
possible  repulsive  interactions.  (33)  is  a  completely  general  way  of  specifying 
any  set  of  2-body  collisions,  and  it  necessarily  describes  invertible  lattice-gas 
dynamics  because  of  the  unitarity  of  the  evolution  operator. 

As  the  dynamics  is  reversible,  the  system  quickly  moves  to  a  maximal  entropy 
state  where  the  net  attractive  interparticle  potential  vanishes.  The  liquid-gas 
coexistence  phase  may  persist  indefinitely  given  a  net  attractive  interaction.  In 
a  reversible  system  a  net  attraction  exists  for  a  short  while,  only  so  long  as 
most  heat  bath  states  are  not  populated.  Once  the  heat  bath  gains  a  significant 
population,  only  the  local  interactions  remain  and  consequent  diffusion  drives 
the  system  back  to  a  disordered  phase.  The  maximal  entropy  state  of  the 
heat  bath  occurs  at  half-filling,  h  =  |,  and  consequently  at  this  heat  bath 
density  it  cannot  encode  any  more  information  about  heating  from  the  lattice- 
gas  so  the  effect  of  the  long-range  interaction  must  become  non-existent.  This 
is  consistent  with  (11),  since  V(d,  h)  =  0  for  h  =  The  heat  bath  population 
exponentially  approaches  its  maximal  entropy  state,  see  figure  7,  starting  from 
a  density  initially  zero;  h(t)  =  |(1  —  e-t/,T)  with  the  observed  time  constant, 
r  =  16.5,  obtained  by  fitting.  The  time  constant,  r,  can  be  increased  by  raising 
the  number  of  heat  bath  states. 
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Figure  7:  Transient  behavior  of  a  reversible  lattice-gas  with  long-range  interactions.  Particle 
density  at  approximately  | -filling.  An  exponential  increase  to  half- filling  with  a  time  constant 
of  16. 5t  is  observed. 
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